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Abstract 

We present a finite difference method to compute the principal eigenvalue and the 
corresponding eigenfunction for a large class of second order elliptic operators including 
notably linear operators in nondivergence form and fully nonlinear operators. 

The principal eigenvalue is computed by solving a finite-dimensional nonlinear min- 
max optimization problem. We prove the convergence of the method and we discuss 
its implementation. Some examples where the exact solution is explicitly known show 
the effectiveness of the method. 
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1 Introduction 

Consider the elliptic self-adjoint operator 

Lu{x) = di {aij{x)dju{x )), (1-1) 

where aij = aji are smooth functions in O, a smooth bounded open subset of M”, satisfying 
> agp for some a > 0. It is well-known that the minimum value Ai in the 
Rayleigh-Ritz variational formula 

. „ - (p(x) Lip(x) dx . J^aij{x)djcp{x)diip{x) dx 

Ai ini Tj jTo mi Tj jT^ 

is attained at some function wi satisfying 

J Lwi{x) + XiWi{x) = 0 X £ Q, 

\ r(;i(x) = 0 X £ dO,. 
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The number Ai is usually referred to as the principal eigenvalue of L in n and wi is 


the corresponding principal eigenfunction. For operators of the form (1.1) and also more 


general linear operator in divergence form there is a vast literature on computational 
methods for the principal eigenvalue, see for example 

General non-divergence type elliptic operators, namely 


Lu{x) = aij{x)diju{x) + bi{x)diu{x) + c{x)u 


( 1 . 2 ) 


are not self-adjoint and the spectral theory is then much more involved: in particular, the 
Rayleigh-Ritz variational formula is not available anymore. In the seminal paper |12j by 
M.D. Donsker and S.R.S. Varadhan, a min-max formula for the principal eigenvalue of a 


class of elliptic operators L including (1.2) was proved, namely 


Ai = - 


Lifix) 
sup— 
</jec2(n),(p>o x&n ^{.x) 


inf 


(1.3) 


In that papers other representation formulas for Ai were also proposed in terms of large 
deviations and of the average long run time behavior of the positive semigroup generated 
by L. A further crucial step in that direction is the paper [6] by H. Berestycki, L. Nirenberg 


and S.R.S. Varadhan, where the validity of formula (1.3) is proved under mild smoothness 


assumptions (H a bounded open set and G (7*^(1!), 6j, c G L°°(n)). Moreover it is 


proved that (1.3) is equivalent to 


Ai := sup{A G M : 3 y? > 0 such that Lip -|- Ay? < 0 in 11}. 

Following this path of ideas, notions of principal eigenvalue for fully nonlinear uniformly 
elliptic operators of the form 

F[u] = F{x, u{x), Du{x), D‘^u{x)) 

have been introduced and analyzed in m, 0 , i, m, m, m- A by now established 
definition of principal eigenvalue is given by 


Ai := sup{A G M : 3 y? > 0 such that F[p] -|- Ay? < 0 in 11} 


(1.4) 


where the inequality in (1.4) is intended in viscosity sense. It is possible to prove under 
appropriate assumptions, see (2.1 )-(|2.2l), that there exists a viscosity solution wi of 


F[rci] -|- Xiwi{x) = 0 
wi{x) = 0 


X G H, 

X G 5H. 


(1.5) 


Moreover the characterization (1.3) still holds in this nonlinear setting. 


As it is well-known, the principal eigenvalue plays a key role in several respects, 
both in the existence theory and in the qualitative analysis of elliptic partial differential 
equations as well in applications to large deviations bifurcation issues [50] , ergodic 

and long run average cost problems in stochastic control |1|. For linear non self-adjoint 
operators and, a fortiori, for nonlinear ones the principal eigenvalue can be explicitly 
computed only in very special cases, see e.g. mm, hence the importance to devise 
numerical algorithms for the problem. But, apart some specific case (see [7j for the p- 
Laplace operator), approximation schemes and computational methods are not available 
in the literature, at least at our present knowledge. 
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The aim of this paper is to define a numerical scheme for the principal eigenvalue 
of nonlinear uniformly elliptic operators via a finite difference approximation of formula 
(1.3). More precisely, denoting by = hZ*^ the orthogonal lattice in M” where h > 0 is 
a discretization parameter, we consider a discrete operator acting on functions defined 
on a discrete subset C Z^ of fl and the corresponding approximated version of (1.3), 
namely 

FhW\{x) 


\h = - inf sup 


( 1 . 6 ) 


As for the approximating operators we consider a specihc class of finite difference 
schemes introduced in HU, HE] since they satisfy some useful properties for the convergence 
analysis. 

We prove that if F is uniformly elliptic and satisfies in addition some quite natural 
further conditions, then it is possible to define a finite difference scheme F^ such that the 
discrete principal eigenvalues Xi^h and the associated discrete eigenfunctions wi^h converge 
uniformly in H, as the mesh step h is sent to 0, respectively to the principal eigenvalue 
Ai and to the corresponding eigenfunction wi for the original problem (1.5). It is worth 


pointing out that the proof of our main convergence result. Theorem |3.2[ cannot rely on 
standard stability results for fully nonlinear partial differential equations, see [3], since the 


limit problem does not satisfy a comparison principle (see Remark 3.1 for details). 


We mention that our approach is partially inspired by the paper m where a similar 
approximation scheme is proposed for the computation of effective Hamiltonians occurring 
in the homogenization of Hamilton-Jacobi equations which can be characterized by a 


formula somewhat similar to (1.3). 


In Section we introduce the main assumptions and we investigate some issues 
related to the Maximum Principle for discrete operators. In Section we study the 
approximation method for a class of finite difference schemes and we prove the convergence 
of the scheme. In Section [4 we show that under some additional structural assumptions 
on Fh the inf-sup problem (1.6) can be transformed into a convex optimization problem 
on the nodes of the grid and we discuss its implementation. A few tests which show the 
efficiency of our method on some simple examples are reported in Section as well. 


2 The Maximum Principle for discrete operators 

We start by hxing some notations and the assumptions on the operator F. Set T = 
n X M X X S'” , where S” denotes the linear space of real, symmetric n x n matrices. 
The function F{x, z,p, r) is assumed to be continuous on T and locally uniformly Lipschitz 
continuous with respect to z, p, r for each fixed x G fl. We will also suppose that the partial 
derivatives Fr, Fp, Fz satisfy the following structure conditions: 

0 < al < Fr < AI, |Fp| < pi, —fjQ < Fz < 0. (2-1) 

for some constants a, A, /tq, pi- A further condition is the positive homogeneity of degree 
1 , that is 

F{x,tz,tp,tr) = tF{x,z,p,r) Vt > 0. (2.2) 

The principal eigenvalue of problem ( |1.5[ ) is defined by 

Al = sup{A : 3 (/? > 0 such that F[ip\ -(- A(/? < 0 in H}, 
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where the differential inequality F\(p\ + Xip < 0 is meant in the viscosity sense. Under 
assumptions (2.1)-(2.2), there exists a viscosity solution of (1.5) and the characterization 
(1.3) of Ai holds (see [8], [TT]1. 

Remark 2.1 It is possible to define 


= sup{A : 3 y? < 0 such that F[ip] + Xip >0 in 11}. 

When F is not odd in its dependence on the Hessian, then in general Ai X^. Of course it 
is possible to see Xf as Ai of some other operator. Hence we will only consider in this paper 
Ai. For example, for the extremal Pucci operators A4'^y^{D‘^u) := supg^j^^^j^jtr{AD‘^u) 
and := mfai<B<Ai H{AD'^u), since M) = the following 

holds 

Remark 2.2 The assumption F^ < 0, i.e. the monotonicity of the differential operator 
in the zero-order term, could be removed. Indeed F := F — cqz, with cq large, satisfies this 
assumption, moreover F and F have the same principal eigenfunction and the eigenvalues 
differ by cq. 


We now describe the discrete setting that we shall consider. Given h > 0, let = /iZ"" 
denote the orthogonal lattice in M”. Let Ff, be a discrete operator acting on functions 


defined in Qh U Z^. We shall consider an approximation of (1.5) (which can be seen also 
as an eigenvalue problem for the discrete operator F^). We look for a number A and a 
positive function w such that 


Fh{x,w{x),[w],c) + = 0 x€llh, 

w{x) = 0 x £ dUh, 


(2.3) 


where 


- h > 0 is the discretization parameter {h is meant to tend to 0), 


- x £ is the point where (1.5) is approximated, 


- w is a real valued mesh function in Z)( meant to approximate the viscosity solution 
of 


- [-ja; represents the stencil of the scheme, i.e. the points in n/i\{x} where the value 
of u is computed for writing the scheme at the point x (we assume that [tcja, is 
independent of w{y) for \x — y\ > Mh for some fixed M £ N). 


We denote by Ch the space of the mesh functions defined on Hh and we introduce some 
basic assumptions for the scheme F^ (see HZ], [HI). 

(i) The operator F^ is of positive type, i.e. for all x £ Ilh, z,t £ M., u,r] £ Ch satisfying 
0 < r]{y) < T for each y £ Pth, then 


Fh{x, z, [u + rffi) > Fh{x, z, [u]^) > Fh{x, z + t,[u + rj],,) 


(ii) 


The operator Fh is positively homogeneous, i.e. for all x £ Clh, z £ M., u £ Ch and 
t > 0, then 

Fh{x,tz, [tu]x) = tFh{x,z, Ma,). 
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(iii) The family of operators {Fh,0 < h < ho}, where ho is a positive constant, is consis¬ 
tent with the operator F on the domain O C M”, i.e. for each u G C'^(n) 

sup \F{x, u{x),Du{x), D‘^u{x)) — Fh{x, u{x), [tt]x)| —^0 as /i —)• 0, 


uniformly on compact subset of Q. 

We study below some properties related to the maximum principle and a comparison 
result for the operator Fh- Let us start by the following definitions: 

Definition 2.1 A function u G Ch is a subsolution (respectively v G Ch is a supersolution) 
of 

Fh{x,u{x),[u]x) = fix) xGith (2.4) 

if 

Fhix, u{x), [u]^) > fix), x gQh 
irespectively, F/j(x,u(x), [v]^) < fix), x G ^h) ■ 

Definition 2.2 The Maximum Principle holds for the operator F^ in fl/j if 


Fhix, uix), [rt]a,) >0 in VLu, 

M < 0 on dVlh, 


(2.5) 


implies rt < 0 in Plh- 


Proposition 2.1 Assume that Fh is of positive type and positive homogeneous and satis¬ 
fies either 


or 


for all z G M, u,rj GCh satisfying 0 < rjiy) and maXy^[.]^r]iy) > 0, 
then Fhix, z, [u + r]]x) > Fhix, z, [u\a:) 


( 2 . 6 ) 


for all z, r G M, u,rj G Ch satisfying 0 < yiy) < r for each y, then 
Fhix, z, [u]x) > Fhix, z-\-T,[u + r]]x) + cqt 


(2.7) 


for some positive constants cq. Then the Maximum Principle holds for the operator Fh in 
Qh ■ 


and M ;= maxj^^ u > Q. Let 
X G be such that uix) = M. Since u < 0 on dPlh-, it is not restrictive to assume that 
there exists y GVlh such that rt(y) < u(x) = M. Hence 


Proof Assume by contradiction that u satisfies (2.5 


0 < Fhix, uix), [n]j;) < Fhix, uix) - M,[u- M]x) 
< Fhix,0, [0]s) = 0, 


a contradiction. A similar proof can be done with the assumption (2.7). 


□ 
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Remark 2.3 The assumptions (2.6) and (2.7) correspond to the uniform ellipticity and, 
respectively, to the strict monotonicity of the operator F with respect to the zero-order 
term. 


The following proposition shows that, as it is known in the continuous case (see for exam¬ 
ple mm, the validity of the Maximum Principle for subsolutions of the operator Fh is 
equivalent to the positivity of the principal eigenvalue for Fh. 

Proposition 2.2 Assume that the scheme Fh is of positive type and that it is positively 
homogeneous. Suppose that for A G M, there exists a nonnegative grid function ip with 
p > 0 in Q,h such that Fh[p] -|- A(/J < 0. If, for t < X, the function u satisfies 

( Fh{x, u{x), [u]x) -\-Tu>0 inHh 

( u < 0 on dflh, 

then u < 0 inTlh, be. = F/i[-] -|- r- satisfies the Maximum Principle. 

Proof Suppose by contradiction that max^^ju} > 0. Let p as in the statement and set 
L{'y) = maxo^ju — 7 (/?} (note that the maximum is taken only with respect to the internal 
points). Then L : [0,oo) —)> M is continuous, decreasing, L(0) > 0 and ^( 7 ) —>■ —00 for 
7 —>■ -|-oo. Hence there exists 7 ' > 0 such that L{'y') = 0. Moreover, since u — 'y'p < 0 on 
dll, we also have max^^ju — 7 V} = 0. Let 0 < 7 < 7 ' be such that 

^A > r (2.8) 

7 

and set = ^p. Then Fh{'fi\ + Xip < 0 and M = maxj^^ju — f/;} = (tt — ’fi){x) > 0 for some 
X G H/j. Hence fi{x) + M = u{x) and tpix) + M > u{x). Since Fh is of positive type, it 
follows that 


Then 


Fh{x, fi{x), [fi]x) > Fh{{x, if{x) + + M]x) = Fh{x, u{x), [fi + M]^ 

> Fh{x,u{x), [u\x). 


Tu{x) > —Fh[u]{x) > —Fh['if']{x) > Xfii{x) = Xjp{x) > X—u{x) 
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and therefore a contradiction to ( 2 . 8 ). 


□ 


The following result gives a comparison principle for (2.4). 


Proposition 2.3 Assume that Fh is of positive type and it satisfies either (2.6) or (2.7). 
Let u and v be a subsolution and respectively a supersolution of (2.4) such that u < v on 
dPth- Then u < v in Llh- 


Proof Suppose by contradiction that M := max^^ju — u} > 0 and let a: G be such 
that u{x) — v{x) = M. Hence v -\- M > u in Qh and it is not restrictive to assume that 
maXj^g[.]_(u + M — u) > 0. It follows that 


f{x) < Fh{x, u{x), [u]x) = Fh{x, v{x) -b M, [u]x) < Fh{x, v{x) -\-M,[v-i- M]^) 

< Fh{x,v{x), [u]s) < f{x) 


and therefore a contradiction. A similar proof can be carried on under assumption (2.7). 

□ 
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3 Approximation of the principal eigenvalue 


In this section we consider a specific class of finite difference schemes introduced in |18j . 
These schemes satisfy certain pointwise estimates which are the discrete analogues of those 
valid for a general class of fully nonlinear, uniformly elliptic equations. 

We assume that for all x G Z^, the stencil [-Jj, of the scheme is given by x + hY where 
Y = {yi,... ,yk} C Z” is a finite set containing all the vectors of the canonical basis of 
M"'. Then we consider a discrete operator in (2.3) given by a finite difference scheme 
written in the form 

Fh[u] = F{x, u, 6hU, 6lu), 


(3.1) 


where F 


X M X 


1 


X JR^ —> M and for y G T, tf G Ch 
u{x ± hy) — u{x) 


h\y\ 


6h,yu{x) = -{S^^yU{x) + 6^yU{x)} = 

Slyuix) = SlydJ^yUix) = 


u{x + hy) — u{x — hy) 
2h\y\ 

u{x + hy) + u{x — hy) — 2u{x) 


/l 2 |y |2 

6hU = {6h,yU : y eY}, 6lu = {Sl^yU : y G T}. 

Set r := M” X M X x and denote by (x, z, q, s) the generic points in T. The operator 
Fh given by (3.1) is of positive type if 


dF 

dsy 

dF 

dz 


dF 

dqy 


>0 yy ^Y, 


< 0 , 


(3.2) 

(3.3) 


and positively homogeneous if 

F{x,tz,tq,ts) = tF{x, z,q, s) Vt > 0. 

Moreover if F in (1.5) satisfies the assumptions then it is always possible to find a 

scheme of type (3.1) which is consistent with F and which, besides (3.2)-(3.3), satisfies for 
all y G T, the bounds 


dF 

l^y| 

dF 

> ao, 

dF 

dF 

dSy 

2 

dqy 

b— - 

dSy 

dqy 


< bo 


(3.4) 


where ao, oq, bo are constants depending on a, A, yo, yi in (2.1) (see [IT], HSj). Note 
that in particular (3.4) implies (2.6). 

We recall some important properties of the previous scheme (for the proof we refer to |18| 1 


Proposition 3.1 Assume (3.2)-(3.4) and let f, g be two given mesh functions. Then for 
every h > 0 sufficiently small there exists a unique solution u/j : —>■ M to the Dirichlet 
problem 


f Fk(x,u(x),[u]x) = f X € Qh, 

[ u = g X G dflh- 


(3.5) 
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Proposition 3.2 Assume (3.2)-(3.4) and let Uh he a subsolution of (3.5). Then 


maxu/j < max 5 + — ■/ h'^\f(x)V 

Q.h ao ^ 

xGilh 


(3.6) 


where the eonstant C is independent of h. Moreover if Uh is a solution of (3.5), then for 
any x,y £Qh 


\uh{x) - Uh{y)\ < C 


k - y\^ 

R 


( 


R 


maxuftH- I h'^\f{x) 

rvcK I < ^ 


Br “0 


(3.7) 


where R = mm{dist(x, SQ/i), dist(x, Sfi/^)}, = B{0,R) n Qh, h and C are positive 

constants independent of h. 


We give an example of a scheme of the form (3.1). Consider the Hamilton-Jacobi-Bellman 
operator 


F{x,u, Du{x), D‘^u{x)) = sup inf L°'^u{x) 


where 


L°‘^u{x) = a“^ {x)DijU + {x)Diu{x) + c°‘^{x)u{x). (3.8) 

It is always possible to rewrite the operator in ( |3.8[ ) in the following form (see m) 

I“^u(x) = af{x)Dl^u -h ¥f{x)Dy^u{x) + c“^(x)?r(x) 


where Dy^u = {Du, yk) and Y = {yi,... ,yk} C Z”" is a finite set containing all the vectors 
of the canonical basis in M"'. Moreover the coefficients 6^^ and satisfy the same 
properties of o“^, and c“^. Then we consider 


Fh[u]{x) := sup inf L?^u{x) 


where 


L'fuix) = {x)5ly^u{x) + bf {x)6h,y,,u{x) + c"^(x)u(x). 

For X G M with H = {1} the previous scheme reads as 

■ r f 0 -/ 9 / su(x + h) + u(x - h) - 2u(x) .u(x + h)-u(x-h) 

sup mf |a /^(x)--+ 6 /^(x)- - - 


(3.9) 

(3.10) 


+ 


+ c' 


,0/3 


(x)u(x)| 


= 0 . 


3.1 The linear case 

In this part we assume that the operator F in is linear, i.e. F[u] = Lu with 

Lu = aij{x)DijU + bi{x)Diu{x) + c(x)u(x) 

and we consider a scheme defined as in ( |3.9| )"( 3T0) ), obviously without the dependence on 
a, j3. 
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Proposition 3.3 Under the assumption (3.4) the eigenvalue problem (2.3) has a simple 
eigenvalue Xi^h £ which corresponds to a positive eigenfunction. The other eigenvalues 
correspond to sign changing eigenfunctions. 

Proof Choose ^ > 0 large enough so that c(x) — ^ < 0 and set 

Lh,i{x, t, [u]x) = Lh{x, t, - ft. 

Let K be the positive cone of the nonnegative grid functions in Ch. For a given grid 


function /, by Proposition |3.1| and Proposition 2^ there exists a unique solution u £ Ch 
to 

Lh,iiy,u{y),[u]y) + f = 0 in 0 / 1 , 

M = 0 on dClh- 

Since Ch is a finite dimensional space it follows that T : Ch ^ Ch defined by T/ = u 
is a compact linear operator. Moreover, if / > 0, then by Proposition |2.1| u > 0 and if 
feK\{0},u = Tf>0. 

Therefore, by the Krein-Rutman theorem [19], r{T) the spectral radius of T is a simple 
real eigenvalue r(T) > 0 with a positive eigenfunction u such that Tu = r(T)u. Hence for 
Ai^/i = r{T)~^ — f, wi = Tu satisfies 

Lh{x,wi{x),[wi]^) + Xi^hWi = 0 in O/i, 

rci = 0 on dXlh. 


□ 


The following characterization of Ai^/i is a simple consequence of Proposition 2.2 


Proposition 3.4 We have 


Xi^h = sup {A : 3 (/9 > 0 s.t. Lh[(p] + Ay? < 0 m H} , 


or, equivalently, 


Xi h = — inf sup 


Lh[T]{x) ] 
t{x) J ■ 


(3.11) 

(3.12) 


Proof Denote by A the right hand side of (3.11). Clearly Ai^/i < A. If Xi^h < ^ then 
there exist /r G (Ai^/i, A) and y? > 0 such that Lh[(p] + < 0. A contradiction follows 

immediately by Proposition 2.2 since the eigenfunction corresponding to Xi^h is positive. 
Hence we have (3.11). 

Let y? > 0 such that L/i[y?](x) + Xip{x) < 0 for x G Clh- Hence 


Consequently 


A < inf 



= — sup 



Ai,/i = sup 
</ 3>0 


— sup 

xaVlh 


ip{x) J j 


= — inf sup 


Lh[T\{x) \ 
t{x) j ■ 


□ 


We give next an upper bound for Xi^h (compare with the corresponding estimate for Ai in 
|6|, Lemma 1.1). 
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Lemma 3.1 Let n = 1 and assume that Br = {|x| < i?} lies in Ll with R < 1. Then 




w 


Proof Given the linear operator 

Lu = a{x)u" + b{x)u' {x) + c(x)m(x), 

let 70, To, h be positive constants such that 70 < a{x) < Pq and \h{x)\, |c(x)| < 6 in 0. 
Let r = R /2 and assume for simplicity that r = Nh for some iV G N. Set = {|x| < r} 
and consider the grid function 

ai = {r'^ - \ih\^f i =-N+ 1 ,... N-I 

Then for i = —N + 1 ,..., — 1 we have 

2 h^^~^ ~ — 4 /ii(r^ — |ihp) + 4 /i^i 

i-2u, ^ 2 _ I ^ 2h2(2z2 + 1). 

h^ 

Denote by a*, and q the coefficients of the linear operator computed at the point x = ih. 
Since h? jir^ — \ih\‘^) < 1 it follows that 

Lh[a]{ih) at _,, , 2 r ^ 

4 crj “ (r^ — \ih\‘^) (r^ — |i/iP )2 * — \ih\‘^) 4 

Pq + hr _ 7ol^^P , b 

~ (r^ — \ih\‘^) (r^ — |i/iP)2 4 


If 


N/iP(7o + Pq + br) > r2(ro + br), 


then the second term in (3.13) dominates the first one and therefore 


In the remaining part of Br, 

Lh[a]{ih) ^ To + br 


4 ( 7 ,: 


LhWKih) ^ b 
4(Ti “ 4 ■ 


+ 7 < 7 H-^(Po + br){-fo + Po + br). 


(r^ — |i/iP) 4 4 7or2 


(3.14) 


(3.15) 


By ( |3.14[ ) and ( |3.15[ ), we get 

Lh[a]{ih)\ ^ C 


sup 

Br 


(Ji 


for i =-N + 1,... ,N - 1. 


To conclude the proof, we show that if for some positive function tp and A G M, Lh[(p\+Xp < 
0, then A < sup 5 ^(—L/i[(T]/fT). For this purpose, assume that A > sup 5 ^(—L/j[(t]/( 7) := r; 
then Lh[a\ +T(7 > 0 in and cr = 0 on dBr, while +t</? < 0. Hence by Proposition 
|2.2[ it follows (7 < 0 in H,., a contradiction, and therefore A < r. □ 
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3.2 The nonlinear case 


We consider now a general discrete operator given by (3.1) and we study the corre¬ 


sponding eigenvalue problem (2.3). In analogy with formula (3.11), we define 


Xi^h = sup {A : 3 > 0 such that Fh[(p] + Ay? < 0} (3.16) 

We prove for each h the existence of a pair {Xi^h,wi,h) satisfying ( |2.3[ ) with wi^h > 0 in 
Qh- 


Proposition 3.5 Assume that Fh satisfies (3.4), / < 0 and X < Xi^h- Then there exists 
a nonnegative solution to 


Fh{x,u{x),[u]x) + Xu{x) = f{x) x£Qh, 

u{x) = 0 X G dVlh. 


(3.17) 


Proof We can assume A > 0, since for A < 0, Fh[u] + Xu satisfies (2.7) and therefore 


by Propositions 3.1 and (2.3) there exists a unique solution to problem (3.17) 


Let us define by induction a sequence Un by setting ui = 0 and, for n > 1 we consider the 
equation: 

Fh{x,Un+l{x), [Un+l]x) = /(x) - XUn, X G 0/,, 

M„+l(x)=0 X G dQ.h. 


(3.18) 


For any n G N there exists a non negative solution Un+i to (3.18). For n = 1, existence 


follows by Proposition [ 3 ^ Moreover since tti = 0 is a subsolution to ( |3.18| ), by Proposition 
2.3 we get U 2 > 0. The existence of a non negative solution at the (n + l)-step is proved 
in a similar way; moreover the solution is non negative since / — Xun < 0. 

We claim now that, for any n > 1, < Un+i- For n = 1 the claim is trivially true since 

U 2 > 0. Assume then by induction that Un > Un-i- Since /(x) — Xun < /(x) — Xun-i it 


follows that Un is a subsolution of (3.18). By Proposition 2.3, we get that Un < Un+i- 


Let us show now that the sequence Un is bounded. Assume by contradiction that it is 
false and set Tin = u„,/|un|oo- Then, by positive homogeneity. Tin is a solution of 


F/j(x,(x), ['iXfi+ijx) — 


fix) 

l^n+1 |c 


-A 


Ur, 


|‘*^n+l|c 


X G n 


h’ 


Since the sequence Un is bounded, then up to a subsequence it converges to a function u, 
while Un/\un+i\oo converges to ku where k = lim„_^oo |^tn|oo/|^tn+i|oo < 1- Hence u > 0, 
|u|oo = 1, = 0 on dQh and 




Fh{x, u{x), [ttja;) + kXu = 0, X GTly^ 

Since 0 < /cA < A and using the fact that for A < Xi^h there exists by definition y? > 0 such 


that Fh[(p] + Ay? > 0 in fl/i, we get a contradiction to Proposition 2.2, Hence the sequence 
Un is bounded, and being in addition monotone, it converges pointwise to a function u 
which solves (3.17). □ 


The next result shows that Xi^h is indeed an eigenvalue for the approximated operator Fh. 
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Theorem 3.1 Assume that satisfies (3.4). Then there exists wi^h > 0 inQ^ satisfying 

(3.19) 


Fhix, wi^h{x), [wi^h]x) + Xi,hWi^h{x) = 0 
wi^h = 0 


X G 

X G dQ.h- 


Moreover that the characterization (3.12) is still valid for the nonlinear operator F^. 


Proof Let Xn be an increasing sequence converging to By Proposition 3.5 there 
exists a positive solution Un of 


Fh{x,Un{x), [Ur, 
Unix) = 0 


+ XnUn — 1 ; 


X G ^h, 

X G dQ.h- 


We claim that Un is not bounded. Assume by contradiction that Un is bounded so that, 
up to a subsequence, Un converges to a function u > 0 which solves 

Fhix,uix),[u\x) + Xi^hU =-I, xG^lh, 
u{x) = 0 X G d^h- 

Then, for e > 0 small enough, u satisfies 

Fhix, uix), [u]x) + iXi,h + £)u = -I + eu < 0 


which gives a contradiction to the definition (3.16). Hence |u„ 
Define now Wn = Un/\un\oo that solves 


oo. 


Fh (x, Wn (x), ['fUn]a:) 4“ XnWn — |nn|oo ’ X G D/j, 

Wnix) = 0 X G dVLh- 

Then, up to a subsequence, Wn converges to a bounded function wi^h which has norm 1 


and which satisfies (3.19), so that wi^h > li¬ 
lt is immediate that (3.12) is still valid for Fh- 


□ 


Remark 3.1 There is a huge literature about the approximation of viscosity solutions 
of first and second order PDFs. In this framework a well established technique to prove 
the convergence of a numerical scheme is the Barles-Souganidis ’method m: besides some 
natural properties of the scheme (stability, consistency, monotonicity), a key ingredient for 
this technique is a strong comparison result for the continuous problem which allows to 
show that a subsolution is always lower than or equal to a supersolution. The comparison 
principle implies in particular that there is at most one viscosity solution of the problem. 


But it is immediate that (1.5) cannot satisfy a comparison principle since w = 0 and the 


principal eigenfunction wi are two distinct solutions of the problem, hence the convergence 
proof cannot rely on the Barles-Souganidis’method and it needs a different argument. 

We now discuss the convergence of the discrete principal eigenvalue to the 
continuous one defined by (1.3). We recall the definition of weak limits in viscosity sense 
(see [3!) 

limsup*u/i(x) := lim sup{u 5 (y) : \x — y\ < h, 6 < h}, 

/i-s-O /1-S-0+ 

liminhu/i(x) := lim inf{M 5 (y) : \x — y\ < h, 6 < h}. 

/i-s -0 /1-S.0+ 
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Theorem 3.2 Assume (2.1)-(2.2), (3.2)-(3.4) and that is consistent with F. Let 
the sequence of the discrete eigenvalues and of the corresponding eigen¬ 
functions, solutions of (2.3). Then ^ —>■ Ai and wi^h wi uniformly in Q as h ^ 0, 

where Ai and wi are respectively the principal eigenvalue and a corresponding eigenfunc¬ 
tion associated to F. 


Proof 

By the positive homogeneity of the scheme, it is not restrictive to assume that maxo^lrci^/i} = 
1, hence the sequence wi^h is bounded. We first prove that 

liminf Ai > Ai. (3.20) 

h—^0 ’ 

Assume by contradiction that liminf/i_,.o Ai^^ = r for some some r < Ai. Consider a 

subsequence, still denoted by Xi^h, such that lim/i_,.o Ai^h = r. Set w = limsup*tci^/i. By 

h^O 

standard stability results in viscosity solution theory, see [3], w satisfies in viscosity sense 

F[t(;] + Tw >0 in (3.21) 


and 


maxtc = 1. 


(3.22) 


Let r/ > 0 be such that for h sufficiently small, Xi^h < t + rj. Hence 
Fh[wi^h] = >-T-r], X G flh- 

Set / = —r — rj, g = 0 and let Uh be the corresponding solution of (3.5), while u is the 
solution of 

J F{x,u{x),Du,D‘^u) = —T — g X € fi, 

\ n = 0 X G dQ. 

Then by Proposition 2.3 and the consistency of the scheme for h sufficiently small 

0 < wi^/i < Ufi < u -G o(l) in (3.23) 


and therefore 


w = 0 on dfl. 


(3.24) 


By (3.21), (3.22) and (3.24) we get a contradiction to the maximum principle for the 
operator F (see ID, mi) and therefore ( |3.20[ ). 

We now prove that 

limsupAi^/i < Ai. (3.25) 

h-s-O 

Assume by contradiction that there exists g > 0 such that 


A := limsup Xi^h > Ai + r/. 

h-5-O 

We consider a subsequence, still denoted by Xi^h, such that lim/i_,.o Ai^/i = A and we set 

w = liminf^tci h- By standard stability results w satishes 0 < w < 1 and 
h^O ’ 


( F[w} + (Ai + g)w < 0 
( w = 0 
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in Ll, 
on dLl, 


(3.26) 













in viscosity sense. Let Xh G be a sequence such that Xh ^ xq & Q and wi^hi^h) = 1 


for all h > 0. By (3.23), xq G 12. We claim that 


w{xo) > 0. 


(3.27) 


Assume by contradiction that uj(xo) = 0, hence there exists a sequence j/h —>■ xq such that 


limh^o m,h(l/h) = 0. By ([3^ with ut = wi^t and / = -Xi^hWyh we get 

/ 

Bn 


\wi^h{xh) - wi^h{yh)\ < C 


\xh - yh\^ 


< C 


R 

Xh - yh\^ 

R 


R 


max wih - \ IA 

R 7 -. ’ ri/f\ I / ^ 


V 


ao 


l,hWlh 





V «0 / 


Since limh^o{wi^h{xh) — wi,h{yh)) = 1 we get a contradiction for h sufficiently small and 


therefore (3.27). 


We are in a position to apply the maximum principle for the continuous problem (see El), 


the definition of Ai. By (3.20) and (3.25) we get lim/i_^o Ai,/i = Ai. 


and so we obtain that ^ > 0. But (3.26) and the positivity of w give a contradiction to 


By (3.7) and a local boundary estimate for wi^h-, see m Thm. 5.1] and [THl Thm.3.], 
we get the equi-continuity of the family and therefore the uniform convergence, 

up to a subsequence, of wi^h to wi with ||tci||oo = 1- The simplicity of the eigenfunction 
associated to the principal eigenvalue Ai gives the uniform convergence of all the sequence 
wi^h to rci. □ 


4 An algorithm for computing the principal eigenvalue 


In this section we discuss an algorithm for the computation of the principal eigenvalue 
based on the inf-sup formula (3.12). In fact we show that this formula results in a finite 
dimensional nonlinear optimization problem. 


4.1 Discretization in one dimension. 


We first present the scheme in one dimension. Note that since the eigenfunction corre¬ 
sponding to the principal eigenvalue vanishes on the boundary of 12/^ and it is strictly 
positive inside, then the minimization in (3.12) can be restricted to the internal points. 
By the formula (3.1) and the homogeneity of J-, we have 


Fh[u]{xi) _ j: ( u{xi + h) - u{xi - h) u{xi + h) + u{xi + h) _ 2 A 
u{xi) V *’ ’ 2hu{xi) ’ h'^u{xi) h‘^) 


We identify the function u{x) with the values t/j, 2 = 0,..., + 1, at the points of the 

grid (with Uq = = 0). Assume that R{x, z,q, s) is linear or more generally convex 

in {q, s). Then the functions G : —>■ defined by 


Gi{x, Ui,..., UnJ = R Ixi, 1, 


Ui—\ UiJ^\ -\- Ui—i 


2hUi 


hm 


/l 2 
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for i = is either linear or respectively convex in f7j+i, Ui-i. Moreover, since 

Ui > 0, G is also convex in Ui. Taking the maximum of the functions Gi over the internal 
nodes of the grid gives a convex function Q : —)> M defined by 

g{Ui,... ,UnJ = max Gi{xi,Ui,... ,UNh) (4.1) 

1=1,...,Nh, 

Hence the computation of Ai^/i is equivalent to the minimization of the convex function Q 
of Nfi variables: this problem can be solved by means of standard algorithms in convex 
optimization. Note also that the minimum is unique and the map is sparse, in the sense 
that the value of Q at Ui depends only on the values at Ui-i and U+i. 

In general, if T{x, z,q, s) is not convex, the computation of the principal eigenvalue is 
equivalent to the solution of a min-max problem in . 

To solve min-max problem we use the routine fminmax available in the Optimization 
Toolbox of MATLAB. This routine is implemented on a laptop and therefore the number 
of variables is modest. A better implementation of the minimization procedure which takes 
advantage of the sparse structure of the problem would allow to solve larger problems. 

Example 1. To validate the algorithm we begin by studying the eigenvalue problem: 

f U)'' + Xw = 0 a: G (0,1), 

\ w{x) =0 a: = 0, 1. 

In this case the eigenvalue and the corresponding eigenfunction are given by 

Ai = TT^, wi{x) = sin(7rx). 


Note that since the eigenfunctions are defined up to multiplicative constant, we normalize 
the value by taking ||rci||oo = ||i(^i,/i||oo = 1 (the constraint for wi^h is included in the 
routine fminmax). Given a discretization step h and the corresponding grid points Xi = ih, 
z = 0,..., Nh + 1, the minimization problem (4.1) is 


Xih = — 


max 

i=l,...,Nh 


Ui-i-i + Ui-i — 2Ui 


hWi 


(with Uq = lAvh+i = 0)- In Table 4.1, we compare the exact solution with the approximate 


one obtained by the scheme (4.1). We report the approximation error for Ai and wi (in 


L°°-norm and L^-norm) and the order of convergence for Ai. We can observe an order of 
convergence close to 2 for Ai and therefore equivalent to one obtained by discretization of 
the Rayleigh quotient via finite elements (see m)- 

Example 2. In this example we consider the eigenvalue problem for a linear 
equation with a discontinuous coefficient 


a{x)w" + Xw = 0 
w{x) = 0 


X G (0,7r), 
X G {0, tt}. 


(4.2) 


where 


a{x) = 


1 for X G [0, 

2 for X G [|^,7r]. 


and k := 


2 +V2 


2V2 


> 1 . 


Proposition 4.1 The principal eigenvalue Ai associated to problem (4.2) is given by k'^ = 
3 + 2 V 2 
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h 

Err(Xi) 

Order{Xi) 

Erroo{wi) 

Err2{wi) 

1.00 • 10-1 

8.0908 • 10-^ 


3.3662 • 10-11 

5.7732 - 10-11 

5.00 • 10-^ 

2.0277-10-^ 

1.9964 

1.4786 • 10-11 

3.8119-10-11 

2.50 • 10“^ 

5.0723 • 10-^ 

1.9991 

6.6613-10-11 

1.8731 - 10-11 

1.25 • 10“^ 

1.2683 • 10“^ 

1.9998 

1.5543-10-11 

6.2524-10-11 

6.25 • 10“^ 

3.1708 • 10-4 

1.9999 

1.2212 - 10-11 

7.1576-10-11 


Table 4.1: Space step (first column), eigenvalue error (second column), convergence order 
(third column), eigenfunction error in L°° (fourth column), eigenfunction error in (last 
column) 


Proof Let 


( \ - j sin(A:x) for x G [0, |^), 

6sin(^+c) for xG[§,7r]. 


'P2 

We choose b and c such that w{0) = w{tt) = 0 and w continuous in Imposing these 
conditions we get 

—^ + c = 7r, and 6sin(—— + c) = 1, 

V 2 2v 2 


i.e. 

n ^ ^ ~ 

c = 7 r(l -—) = 7 r (---). 

Furthermore, using that + c = |, we get 

TT TT 

lim w'{x) = kcos{—) = 0, lim w\x) = bkcos {—^ + c) = 0, 


2k 


' 2 V 2 


hence w G (^^([O, vr]). 

On the other hand, since w is not in we show that it satisfies (4.2) in the sense of 
viscosity solutions. For any {p,q) G we get p = 0 and q > This implies 

that for both a = 1 and a = 2: 

aq > — 

so tc is a subsolution. For any (p, q) G we get p = 0 and q < —k'^. This implies 

that for both o = 1 and a = 2: 

TT 


aq < -k w{ — ), 


so rc is a supersolution. 


□ 


In Table 4.2, we compare the exact solution with the approximate one obtained by means 
of the scheme 


Xi h = — mill 

’ um^h 


max 

i=l,...,Nh 


a{ih) 


Ui+i + Ui-i — 2Ui 

TO 


(with Uq = = 0). The rates are not very good, but the problem is out of our setting 

since F is discontinuous and the error is very sensible to the chosen grid. In Figure 4.1, 
we report the graph of the exact and approximate eigenfunctions for h = 0.1. 
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h 

Err{Xi) 

Order{Xi) 

Errooiwi) 

Err2{wi) 

0.1571 

0.1197 


0.0213 

0.0563 

0.0785 

0.0476 

1.3303 

0.0090 

0.0383 

0.0393 

0.0347 

0.4576 

0.0065 

0.0391 

0.0196 

0.0157 

1.1417 

0.0030 

0.0264 

0.0098 

0.0061 

1.3596 

0.0012 

0.0149 


Table 4.2: Space step (first column), eigenvalue error (second column), eigenfunction error 
in L°° (fourth column), eigenfunction error in (last column) 



Figure 4.1: 


Exact and approximate eigenfunctions for /i = 10 ^ 


Example 3. The Eucik spectrum of A is the set of pairs (/r, a/r) G for which 
the equation 

—Am = fiu'^ — a^u~ 

has a non-zero solution, where u'^{x) = max{M(x),0} and u~{x) = max{—M(a:),0}. Eor 
fixeda > 0 the previous problem is equivalent to 

mini Am, — Am} + fiu = 0 if a > 1, 
a 

maxi Am, — Am} + fiu = 0 if a < 1. 
a 


17 


















For details see m- Hence the Fucik spectrum can be seen as the spectrum of a nonlin¬ 
ear operator involving the maximum or minimum of two linear operators. To find the 
corresponding principal eigenvalue we apply the scheme ( |4.1[ ). In Table 4.3, we report 
the corresponding approximation error for Ai in the case a = 1/2 and H = [ 0 , 7 r] (by the 
convexity of the solution the eigenvalue for the continuous problem coincides with the one 
of the second derivative operator in [ 0 , vr] i.e. Ai = 1). 


h 

Err{Xi) 

Order{Xi) 

1 . 00 - 10-1 

0.0809 


5.00-10-^ 

0.0203 

1.9964 

2.50-10-^ 

0.0051 

1.9991 

1.25 • 10“^ 

0.0013 

1.9998 

6.25 • 10“^ 

0.0003 

2.0000 


Table 4.3: Space step (first column), eigenvalue error (second column), convergence order 
(third column) for the Fucik spectrum with 0 = 1/2 


Example 4. Consider the eigenvalue problem for the p-Laplace operator 

Dwi) + = 0. 


This example does not fit exactly in the framework of this paper since the operator is not 
uniformly elliptic. However, the following formula 




inf sup 


Fh,p[^](y) ] 
(p{y)P-^ J 


where Fh^p is a finite-difference approximations of Fp produces a good approximation of 
the principal eigenvalue of the p-Laplace operator in the interval (a, b) whose exact value 
is given by 

^ ^ ( 6 -a)psin(^)' 

In Table 4.4 we report the approximation error and the corresponding order of convergence 
for the principal eigenvalue of the p-Laplace operator for p = 4 (in this case A 4 ~ 73.0568). 
It is also known (see m) that if H is a ball, the eigenfunction Wp corresponding to the 


h 

Err[X/C) 

Order{Xi) 

1 . 00 - 10-1 

2.6770 


5.00-10-^ 

0.6210 

2.1079 

2.50-10-^ 

0.1457 

2.0912 

1.25 - 10-^ 

0.0347 

2.0724 

6.25 - 10-^ 

0.0083 

2.0581 


Table 4.4; Space step(first column), eigenvalue error (second column), convergence order 
(third column) for the p-Laplace operator with p = 4 

eigenvalue Ap converges for p —)■ 00 to d{x,dXl). In Figure 4.2, we draw approximations 
of Wp computed by the scheme for various values of p and we observe the convergence of 
these functions to d{x, { 0 , 1 }) for p increasing, as expected by the theory. 
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Figure 4.2: Approximate eigenfunction Up^h for p = 2,4, 6,8,10 and /i = 10 ^ 


4.2 Discretization in higher dimension. 

We now consider the eigenvalue problem in Arguing as in the 1-dimensional case we 
write 


u{x) 

for i = 1,. 


= F 


{x,: 


u{x + hy) — u{x — hy) \ 


2h\y\u{x) 


7 


VdY 


f u(x + hy) -\- u(x — hy) 2 '[ ^ 

1 h^\y\Mx) 


(4.3) 

,Ny and Fy defined as in (3.1), Y the stencil and Ny the cardinality of 


^h- Hence if the function F{x, z,{qy}y^Y,{sy}ye:Y) is linear or more generally convex in 
the variables qy and Sy, y ^ Y, then the computation of the principal eigenvalue Xi^y is 
equivalent to the minimization with respect to the vector U G of the convex function 
Q : —)• M obtained by taking the maximum with respect to x G in (4.3). Therefore 

this problem can be solved by means of some standard algorithms in convex optimization. 

Example 5. Consider the problem 


f Aw + Xw = Q xG(0, 1)^, 

\ w(x) = 0 X G d((0,1)'^). 


The eigenvalue and the corresponding eigenfunction are given by 

Ai = 27r^, w(xi,X 2 ) = sin(7rxi) sin(7rx2) 
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(the eigenfunctions are normalized by taking ||tc||oo = 
five-point formula for the discretization of the Laplacian. 


= 1). We use a standard 
In Table 4.5, we compare the 
We report the 


exact solution with the approximate one obtained by the scheme (4.1). 
approximation error for Ai and wi (in L°°-norm and L^-norm) and the order of convergence 
for Ai. We can observe an order of convergence close to 2 for Ai and therefore equivalent 
to one obtained by discretization of the Rayleigh quotient via finite elements, see m- 


h 

Err[\i) 

Order{\i) 

Errooiw) 

Err2{w) 

2.00 • 10“^ 

0.4469 


0.0801 

0.2256 

1.00 • 10“^ 

0.1338 

1.7397 

0.0203 

0.1137 

5.00 • 10-^ 

0.0368 

1.8629 

0.0056 

0.0590 

2.50 • 10“^ 

0.0097 

1.9297 

0.0015 

0.0301 


Table 4.5: Space step (first column), eigenvalue error (second column), convergence order 
(third column), eigenfunction error in L°° (fourth column), eigenfunction error in (last 
column) 

Example 6. We consider the eigenvalue problem for the Ornstein-Uhlenbeck 
operator 

Aw — X ■ Dw + Ate = 0, X G (—1,1)^ 

with homogeneous boundary conditions. The eigenvalue and the corresponding eigenfunc¬ 
tion are given by 

Ai = 4, t(;(xi,X 2 ) = (1 - x?)(l - X 2 ), 

with the eigenfunctions normalized by taking ||rc||oo = ||?i'i,/i||oo = 1- The Laplacian is 
discretized by a five-point formula. In Table 4.6, we report the approximation error for Ai 
and the corresponding order of convergence. 


h 

Err[\i) 

Order(Ai) 

4.00-10-1 

0.1524 


2.00-10-1 

0.0392 

1.9592 

1.00-10-1 

0.0103 

1.9250 

5.00-10-^ 

0.0027 

1.9580 


Table 4.6: Space step (first column), eigenvalue error (second column), convergence order 
(third column) 


References 

[ 1 ] S.N. Armstrong, The Dirichlet problem for the Bellman equation at resonance. J. 
Differential Equations 247 (2009), 931-955. 

[2] I. Babuska, J. E. Osborn, Finite element-Galerkin approximation of the eigenvalues 
and eigenvectors of selfadjoint problems. Math. Comp. 52 (1989), 275-297. 

[3] G. Barles, P. Souganidis, Convergence of approximation schemes for fully nonlinear 
second order equations. Asymptotic Anal. 4 (1991), 271-283. 


20 





















[4] A. Bensoussan, Perturbation methods in optimal control. Wiley/Gauthier-Villars Se¬ 
ries in Modern Applied Mathematics. John Wiley Sz Sons, Ltd., Chichester; Gauthier- 
Villars, Montrouge, 1988. 

[5] H. Berestycki, I. Capuzzo Dolcetta, A. Porretta, L. Rossi, Maximum Principle and 
generalized principal eigenvalue for degenerate elliptic operators. J. Math. Pures Appl. 
103 (2015), no. 5, 1276-1293. 

[6] H. Berestycki, L. Nirenberg, S.R.S. Varadhan, The principal eigenvalue and maximum 
principle for second-order elliptic operators in general domains. Comm. Pure Appl. 
Math. 47 (1994), no. 1, 47-92. 

[7] R.J. Biezuner, G. Ercole, E.D. Martins, Computing the first eigenvalue of the p- 
Laplacian via the inverse power method. J. Eunct. Anal. 257 (2009), no. 1, 243-270. 

[8] I.Birindelli, F. Demengel, First eigenvalue and maximum principle for fully nonlinear 
singular operators. Adv. Differential Equations 11 (2006), no. 1, 91-119. 

[9] I. Birindelli, F. Leoni, Symmetry minimizes the principal eigenvalue; an example for 
the Pucci’s sup operator. Math. Res. Lett. 21 (2014), no. 5, 953 - 967. 

[10] D. Boffi, Finite element approximation of eigenvalue problems. Acta Numer. 19 

( 2010 ), 1 - 120 . 

[11] J. Busca, M. Esteban, A. Quaas, Nonlinear eigenvalues and bifurcation problems for 
Pucci’s operator. Ann. Inst. H. Poincare Anal. Non Lineaire 22 (2005), no. 2, 187-206. 

[12] M.D. Donsker, S.R.S. Varadhan, On the principal eigenvalue of second-order elliptic 
differential operators. Comm. Pure Appl. Math. 29 (1976), 595-621. 

[13] D. Gomes, A. Oberman, Computing the effective Hamiltonian using a variational 
approach. SIAM J. Control Optim. 43 (2004), 792-812. 

[14] W. Huang, Sign-preserving of principal eigenfunctions in PI finite element approxi¬ 
mation of eigenvalue problems of second-order elliptic operators. J. Comput. Phys. 
274 (2014), 230-244. 

[15] H.Ishii, N.Ikoma, Eigenvalue problem for fully nonlinear second-order elliptic PDE 
on balls. Ann. Inst. H. Poincare Anal. Non Lineaire 29 (2012), 783-812. 

[16] P. Juutinen, P. Lindqvist, J.J. Manfredi, The oo-eigenvalue problem. Arch. Ration. 
Mech. Anal. 148 (1999), no. 2, 89-105. 

[17] H. J. Kuo, N. Trudinger, Linear elliptic difference inequalities with random coefficients. 
Math. Comp. 55 (1990), no. 191, 37-53. 

[18] H. J. Kuo, N. Trudinger, Discrete methods for fully nonlinear elliptic equations. SIAM 
J. Numer. Anal. 29 (1992), 123-135. 

[19] M. G. Krein, M. A. Rutman, Linear operators leaving invariant a cone in a Banach 
space. Amer. Math. Soc. Translation (1950), no. 26, 128 pp. 

[20] P.-L. Lions, Bifurcation and optimal stochastic control. Nonlinear Anal. 7 (1983), 
177-207. 


21 



[21] C. Pucci, Maximum and minimum first eigenvalues for a class of elliptic operators. 
Proc. Amer. Math. Soc. 17 (1966), 788-795. 

[22] H. F. Weinberger, Variational methods for eigenvalue approximation. Conference 
Board of the Mathematical Sciences Regional Conference Series in Applied Math¬ 
ematics, no. 15, SIAM, Philadelphia, 1974. 


22 



